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Abstract 

A medium which is an isotropic chiral medium from the perspective of a co-moving observer 
is a Faraday chiral medium (FCM) from the perspective of a non-co-moving observer. The 
Tellegen constitutive relations for this FCM are established. By an extension of the Beltrami 
field concept, these constitutive relations are exploited to show that planewave propagation is 
characterized by four generally independent wavenumbers. This FCM can support negative 
phase velocity at certain translational velocities and with certain wavevectors, even though the 
corresponding isotropic chiral medium does not. The constitutive relations and Beltrami-like 
fields are also used to develop a convenient spectral representation of the dyadic Green functions 
for the FCM. 

Keywords: Beltrami field, Bohren transform, dyadic Green function, Faraday chiral medium, 
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1 Introduction 



A fundamental issue in electromagnetics is the variation in the perceived properties of a linear 
medium according to the observer's inertial frame of reference. Interest in this topic dates from 
the earliest days of the special theory of relativity and it remains an active area of research. 
Electromagnetic fields in mediums which are isotropic dielectric-magnetic from the perspective 
of a co-moving observer have been widely studied (Chen 1983; Kong 1986; Pappas 1965). Recent 
studies involving an isotropic dielectric-magnetic medium have demonstrated that planewave 
propagation with positive, negative or orthogonal phase velocity can arise depending upon the 
observer's inertial frame of reference (Mackay & Lakhtakia 2004a; Mackay et al. 2006). The 
electromagnetics of simply moving plasmas have also been extensively investigated (Chawla &; 
Unz 1971). 

In this paper we consider electromagnetic fields in linear isotropic chiral mediums moving 
at constant velocities. A natural formalism for investigating the electromagnetic properties 
of an isotropic chiral medium, as observed from a co-moving inertial frame of reference, is 
provided by Beltrami fields. The defining characteristic of a Beltrami field is that the curl 
of the field is a scalar multiple of the field itself (Lakhtakia 1994a, b). These fields are useful 
for analysis of a wide range physical phenomenons, as in astrophysics (Chandrasekhar 1956, 
1957), hydrodynamics and magnetohydrodynamics (Dritschel 1991), thermoacoustics (Ceperley 
1992), chaotic flows (McLaughlin and Pironneau 1991), plasma physics (Yoshida 1991) and 
magnetostatics (Marcinkowski 1992), for example. In the following sections, an extension of the 
Beltrami field concept is developed to investigate the electromagnetic properties of an isotropic 
chiral medium as observed from a non-co-moving inertial frame of reference. 

In earlier studies involving isotropic chiral mediums moving at a constant velocity, the 
Lorentz-transformed wavevector and Lorentz-transformed angular frequency have been utilized 
to explore Doppler shift and aberration (Engheta et al. 1989); the scattering response of an elec- 
trically small sphere made of an isotropic chiral medium has been formulated (Lakhtakia et al. 
1991); and planewave propagation has been investigated for relatively low translational speeds 
(Hillion 1993; Ben-Shimol & Censor 1995, 1997). Reflection and transmission coefficients for a 
uniformly moving isotropic chiral slab have also been calculated using the Lorentz-transformed 
electromagnetic fields (Hinders et al. 1991). 

In contrast to these earlier works, the analysis presented in the following sections begins 
with a derivation of the Tellcgen constitutive relations, from the perspective of a non-co- 
moving observer, for a medium which is an isotropic chiral medium for a co-moving observer. 
By means of the Bohren transform and the consequent introduction of Beltrami-like fields, 
these constitutive relations are exploited to consider planewave propagation — specifically, the 
propensity for negative phase velocity — from the perspective of a non-co-moving observer. The 
constitutive relations, combined with Beltrami-like fields, are also used to develop a convenient 
spectral representation of the dyadic Green functions for the isotropic chiral medium moving at 
constant velocity. 

As regards notational matters, 3-vectors are underlined whereas 3x3 dyadics are double 
underlined. The identity 3x3 dyadic is written as /. Vectors with the "symbol overhead are unit 
vectors. The operators Re { • } and Im { • } deliver the real and imaginary parts, respectively, of 
complex-valued quantities; the superscript * denotes a complex conjugate; and i = ^/—T. The 
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permittivity and permeability of free space are denoted eo and /Xo, respectively; Cq = l/y^e^/x^ is 
the speed of light in free space. 



2 Constitutive relations 



We contrast the electromagnetic properties in two different inertial frames of reference, denoted 
as S' and S, where S' moves with constant velocity v = vv relative to S. The spacetime 
coordinates (r', t') in frame S' are related to the spacetime coordinates (r, t) in frame S by the 
Lorentz transformation 



wherein 



y = / + (7- l)vv, 
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(1) 

(2) 



and the relative speed (5 = v/co. 



In the absence of sources, the spatiotemporal variations of the (time-domain) electromagnetic 
fields in frame S' are related as 



V'xl'(r',O + ^^'(r',O = 

v'x^V,0-|7^'(2:',0 = o J 



whereas those in frame S are related as 



d ^ 



(3) 



V X E{r, t) + -^B{r, t) = Q 

V xH{r,t)-^D{r,t) = 



(4) 



as dictated by the Lorentz covariance of the Maxwell curl postulates. The primed and unprimed 
fields in (3) and (4) are connected via the Lorentz transformation as (Chen 1983) 



Mil!, t') = 7 \xr^ ' Ml, t) + vx b{l, t)] , 

• B{r, t) - \ vx E{r, t) 



B:{r',t') 



Y- 



M!{r', t') = 7 \Y-^ ' H{r, t)-vx D{r, t)] , 



M{r',t') = 7 



XT^ • D{r, t) + ^vx H{r, t) 



(5) 

(6) 

(7) 
(8) 



Let us consider a homogeneous medium, which is an isotropic chiral medium from the per- 
spective of an observer co-moving relative to the frame S'. In the most general scenario, the 
medium is spatiotcmporally nonlocal. From the perspective of the observer co-moving rela- 
tive to S', the constitutive relations of the medium may be expressed in the Tellegen form as 
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(Lakhtakia 1994b) 



— ^0 e'{s ,u') (r — s ,t' — u) duds' 

Js' Ju' 

+i^/e;^, [ [ |'(s', u') H^ir! - s', t' - u') du'dd, 

Js! Ju' 

S'(r', t') = -i^/I^o [ [ f U', u') ^'(r' - s', t' - u') du'ds' 

Js' Ju' 

+Ho [ [ /i'U', u') H^ir! - s', t' - u') du'ds!, 



(9) 



(10) 



with the real-valued (time-domain) constitutive parameters e'(r:',t'), $,'{r',t') and jl'ir'jt'). By 
implementing the spatiotemporal Fourier transforms (Lakhtakia Sz Weiglhofer 1996) 



(11) 







' E!ir',t') ' 






iri It' 




yex.p[i{Lo't' -k! •/)] dt'dr! 






. B!{r',t') . 





and 



4(ifc',c^') 




r' Jt' 



X'iz', t') exp [i {u't' - k! ' r')] dt'dr', 



(X = e,e,//), (12) 



along with the convolution theorem (Walker 1988), the frequency-domain constitutive relations 
emerge as 



^J(fc', cv') = eo ej(fc', cv') g0,io') + zV^Ctt(^', ^') E^^', 
B!0, a;') = -zv^ejl^', ^0 :^J(fc', u') + /Xo fi'^{k', u') ^J(fc', u') 



(13) 



In many applications the effects of spatial nonlocality are negligible in comparison to those 
of temporal nonlocality. The constitutive relations (13) may then be approximated as 



DUl'^u') = e,e'{u')E:{r!,uj') + i^,i'{J)H!{r!,u') 
B'{r',u;') = -i,/^oe{^')E'ir',u;') + fi,ti'{u;')H'{r',u') 

for spatially local mediums, wherein 



(14) 



' E!{r',Lo') 
H!{r!,uj') 
D!{r',uj') 

[ B'{r',u') 



' l'(r',t') 



11^;;;) ^exp(..Y)..' 

I S'(r:',t') 



(15) 



and 



i'{t')ex.-p{ioj't') dt', 
with x'i't') = x'd',*') for X = e, ^ and /x. 



(X = e,?,M), 



(16) 
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Let us now proceed to develop the frequency-domain constitutive relations in the reference 
frame S. After using (5), (7) and (8) to substitute for E_ {r^,t'), H_ {r^ ,t') and {r!,t'), respec- 
tively, the constitutive relation (9) may be expressed in terms of S fields as 



D{r,t) 




e'(l', u')E{r - Z • - -fvu', t-u'-f- H-L^] 



+Y • iv X 



e'{s', u')B{r - H • s' - 7W, t-u'-i 



7s • 



du'ds' 



f (s', u')H{r - y • s' - 7W, t-u'^- ' - ) 



-Y • Ivx 



^'{s', u)D{r — H • 1' — jvu , t — u'j 



7s • V 



-Y 



\v X H_{r, t) 

,''0 



Similarly, the constitutive relation (10) may be expressed in terms of S fields as 



duds' 



(17) 



B{r,t) = -^^/^o / \ i' {s! ,u')E{r- Y ' s! - -fvu' ,t - u'j - 



7 s • V, 



+Y ' <v X 



I'd', u')B{r -Y-s'- -fvu, t-u'-f 



7£ * 



duds' 




il'{s', u')H{r - F • 1' - 7W, t-u'^- ' - ) 



— y 'ivx 



fi'W, u')D{r -Y's'- jvu', t-u'j- ' " ) 



\v X E{r, t) 



du'ds' 



(18) 



by using (5), (7) and (6) to substitute for E_ {r!,t'), H_ {r',t') and 5 {r',t'), respectively. Imple- 
mentation of the spatiotemporal Fourier transforms 



r E0,u) ] 




r Mr,t) ] 






mr,t) \ 






Dir,t) 1 









exp [i {ut — * r)] dt dr 



(19) 



and (12) with (17) delivers 



D^{k,^) = eoe'0,u'){E^{k,u;)+Z' [v x R^ik, 00)]} 

+iV^o^i{k',^') {H^{k,u^) - ^ • [vx D^{k,uj)] } 

-Y ' i,vxH^{k,oj) 



(20) 
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and with (18) yields 



+^o/i«(fc', w') {H^{k,io)-Y' [v X D^ik,u;)] } 



+Y 



-^vxE^{k,u) 



(21) 



In the derivation of (20) and (21), the principle of phase invariance (Pappas 1965; Kong 1986) 
has been invoked to obtain the relations 



fe' = y • fc — ^ V, Lo' = J {uj — k ' v) . 



(22) 



The two independent expressions (20) and (21) which relate D^{k, oo) and B^f^{k, u) to E^{k, lo) 
and H.^{k,oj) can be manipulated to deliver the S frequency-domain constitutive relations 



(23) 



D^ik, u) = eog^ik, uj') • E^ik,uj) +i./e^l^^{k,uj') • H^{k,uj) 
E^{k,i^) = ' E^ik,uj) + lJ-o^^{k',uj') ' H^{k,uj) 

Herein the 3x3 constitutive dyadics all have the form 

K^ik',u;') = xl{k',u;')L- ixl{k',io')v x 1+ [xf(fc',c^') - x\{k',io')] vv, (x = e,^,//). (24) 

The unprimed constitutive parameters emerge as 

4(fe',a;') = e'0,u;') {^^ [e'0 ,u;'mk' " (^',a;')] - l} {P' - l) Mk','^') ] 
eUk!,u;') = 2e'0,J)i'0,J)l3 (l - /J^) A«(fc',a;') 



) , (25) 



ej(fc',a;') = Cj(fc',a;') {(3^ ,u'mk' ,u;') - ej'(^',a;')] + 1} (l - P') ^iik',^') ) 
q{k!,io') = ^(^e'0,J)iJi'0,io') - 1] [1 - f5h\{U,J)iJi'0,J)\ 

~i'i\k!,^') - 2e'^{k!,u;')^i'0,u') - l] - 1} j A„(fe', a;') 

q{k',u;')=e^{k:,u;') 

M§(fe',a;') = f,'0,u;') {[3' [e'0 ,u;')fi'0 ,u;') - ^^^(fe',^;')] - l} (Z?' - l) A„(A;',a;') 
M^(fe',a;') = 2^^l{k',u;')ei{k',u;')P (l - P^) A^ik^u') 
fil{k!,co') = li'^{k',u;') 



>,(26) 



,(27) 



with 



1-2/32 [4(fc',u;')/.J(A;',a;')+^f(Ai',u;')] 
[e'0,u'mk:,u') - ^'i\k\uj'f. 
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(28) 



The constitutive relations (23) in S reduce to those in E' specified by (13) in the hmit v ^ 0. 
Wc note the similarly of these constitutive relations (23) to the tensor formulation that is 
used in plasma physics (Melrose & McPhedran, 1991). Interestingly, the constitutive dyadics 
(24) have the same form as that ascribed to a Faraday chiral medium (FCM) (Weiglhofer & 
Lakhtakia 1998). Hitherto, FCMs have been conceptualized as homogenized composite mediums 
(Weiglhofer & Lakhtakia 1998; Engheta et al. 1992), arising from blending together an isotropic 
chiral medium with either a magnetically biased ferrite (Weiglhofer et al. 1998) or a magnetically 
biased plasma (Weiglhofer k. Mackay 2000) . Through the homogenization process, the natural 
optical activity of isotropic chiral mediums (Lakhtakia 1994b) is combined with the Faraday 
rotation exhibited by gyrotropic mediums (Lax & Button 1962). 

For the spatially local medium represented by the S' constitutive relations (14), the corre- 
sponding constitutive relations in E are delivered from (23) as 

D{r, w) = eo |(a;') • E{r, lj) + i^/e^^{Lj') - H{r, u) 
B{r, uj) = -i,/e^^{uj') • E{r, w) + /Xq /x(a;') • H{r, to) 

with the 3x3 constitutive dyadics defined as in (24), but with no dependency on fc'. 

As an illustrative example, let us consider the case of the medium specified by the S' con- 
stitutive parameters e' = 6.5 + il.5, ^' = 1 + iO.2, and ^' = 3.0 + iO.5. The corresponding 
constitutive parameters in E, as specified in (25)-(27), are plotted in Figure 1 against the 
relative speed /3 G [0,1). The parameters are constrained such that Re{x*,x^} R-ejx'}, 
Im|x*;X^} ~^ I™{x'} and \x^\ — > as /? ^ for x = and ^. In Figure 1, whereas 
Re{x*,X^} a-iid Im|x*)X^} become vanishingly small as (i approaches unity for X = ^ 
jjL, as do Re {4*} and Imj^*,^^}, this is not the case for Re {4^}. The parameters x^ a^'e 
independent of j3 for x = A* and ^. 

3 Planewave propagation 

Let us now consider the propagation of plane waves in spatially local mediums of the chosen 
kind. A plane wave characterized in frame S by the wavevector k and angular frequency oj is 
related to a plane wave characterized by the wavevector fc' and angular frequency to' in frame 
E' by the relations (22). The Doppler shift and aberration arising from the transformation 
from S' to E have been explored previously (Engheta et al. 1989). In the remainder of this 
section, we exploit the frequency-domain constitutive relations (14) and (29) for the two frames 
to consider planewave propagation and, in particular, investigate the phenomenon of negative 
phase velocity for the isotropic chiral medium moving at constant velocity. A central element 
in the analysis is the introduction of Beltrami and Beltrami-like fields. 

Planewave propagation in reference frame E', with field phasors of the form 

^'(r',Lj') =^o(w') exp(i^' -r') 
H!{r',u') = H!q{u') exp(ife' • r') 

is a well-documented matter (Lakhtakia 1994b). It is both mathematically expedient and 
physically insightful to implement the Bohren transform and introduce the Beltrami field phasors 



(29) 
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(Bohren 1974) 



Q'i(£',a;') = X [E:{r',u;') + ir,'{u;')H'ir',u;') 



2 



with the intrinsic impedance 



n'ico') 



//o/i'(w') 



.1/2 



eoe'(^') 

Thereby, the frequency-domain Maxwell curl postulates in reference frame S', namely 

V X ^'(r', Lj') - iLo'gir;, Lo')=0 I 

V X H'{r',io') + iu;'D'{r',u;') = Q J ' 

may be recast as two uncoupled first-order differential equations, which yield 



ik! X q; (r', Lo') + i-iYk'xo;^ oj'^ 



(£ = 1,2), 



(31) 



(32) 



(33) 



(34) 



for plane waves (30). Regardless of the direction of propagation, two wavevectors A;' G {^} with 
= k'^k^ and corresponding wavcnumbcrs fc' = k[^k'^ are supported, where (Lakhtakia 1994b) 



k[ = ^e'{u;')fi'{u')-e{u;') 
k'^ = ^e'{u;')f,'{u')+e{co') 



(35) 



and kg = uj' y/e^pLQ. 



The nonreciprocal bianisotropic nature (Krowne 1984) of the medium specified by (29) in 
reference frame S leads to more complicated planewave characteristics than in S'. Following 
the strategy used for frame S', it is helpful to utilize the field phasors 



1 



id,^) = [Mr,uj) + ivl {J) H{r,uj) 



> . 



(36) 



This enables the frequency-domain Maxwell curl postulates in frame S, namely 

V X E{r, uj) - iojB{r, w) = 1 



V X H{r, uj) + iujD{r, w) = 



to be decoupled as 



ikxQ^{r,u:) + {-lfK~^{uj')'Q^{r,uj) = Q, 



(£=1,2), 



for plane waves 



^{z) = E.Q cxp{ik • r) I 
IL{l) = ILo exp(i A; • r) J 
with wavevector k = kk and ko = to^/ei^. The 3x3 dyadics in (38) are given as 

g{uj') = 4{uj')L-iKi{io')vxl+ [k'^{u')-kl{u')]vv, (£=1,2), 



(37) 

(38) 
(39) 

(40) 



with 



6' (a;') J 



1/2 



1/2 



1/2 



e-(a;0+r(^')^ 



LM'(a;')J 
6'(a;') 



/x-(u;')-r(u;0 



1/2 



>, {iy = t,g,z). (41) 



The unprimed constitutive parameters in (41) are defined as in (25)~(27), but with no depen- 
dency on A;'. Whereas Q'^ fc', w') are Beltrami field phasors (Lakhtakia 1994a), (l,'*^) should 
be called Beltrami-like field phasors. 



The dispersion relations 



det 



ik X /+ {-lYk„^{uj') 



0, 



(£=1,2), 



(42) 



arise immediately from (38). For an arbitrary direction of propagation specified by the relative 
orientation angle 9 = cos~^ v • k, the dispersion relations (42) may be expressed as the pair of 
quadratic equations 



aek +bek + ce = 0, 



{£=1,2) 



(43) 



wherein the relative wavenumber k = k/ko- The coefficients in these equations are given as 



bi = 2/3cos^ jl- 



cos^ 9 + sin^ I 



y^V)7R)-C'(c^') 



V^V)7M+C'(^') 



(44) 



and 



02 



-l + /?2( <j y/e'{Lu') ii'{lu') - e'(w')J cos^^ + sin^e 

2 



X 1-/32 v'iv)MV)+e'(c^') + v^v)mV)+4'('^') 



62 = 2/3 cos 61 



Ve'(a;')/x'(a;') + C'(^') 



n2 



-1-1+1 



C2 



(45) 



Thus, four relative wavenumbers k G |^i, ^2, ^3, emerge as the roots of (43). While these 
may be straightforwardly extracted from (43), explicit algebraic representations of the wavenum- 
bers are generally cumbersome. The following two special cases are noteworthy exceptions. For 
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propagation parallel to the direction of translation (i.e., ^ = ^) we have 

~ _ ViVFM (1 /5 [1 - e'{u;')^l'{J) + ei^')] - (1 + 



I -13 {(3 [e'{Lo')^l'{,o') - e'2(^0] - 2e'(c^')} 



ko = 



h = 



^e'{u:')^i'{u:') (l - /g^) + /? [l - e' {lo') + ^W)] + ^^(c^O (l + 
1 - [e'(^')/"'(^') - ^2(^0] - 2e'(c^')} 
V€^(a;OM^(a;0 (l - /g^) + /3 [l - 6^(.;0a^^(.;0 + + ^'juj') (l + /?^) 

1 - [e'(^')/^'(^') - - 2^'(^')} 



A;4 = 



- ViVFR) (1 /3 [1 - + ei^')] - (1 + 



whereas the relative wavenumbers are delivered as 



(46) 



k2 = < 



1 2 



1-/32 



1-^2 



1/2 



1/2 



(47) 



k3 = -ki 
k4 = —k2 



for propagation perpendicular to the direction of translation (i.e., k • v = 0). 

By way of numerical illustration, let us return to the constitutive parameters used for Figure 1. 
The corresponding relative wavenumbers in E, computed as the roots of (43), are plotted in 
Figure 2 against relative speed /3 G [0, 1) and wavevector orientation angle 9 G [0, tt). For clarity, 
the wavenumbers in Figure 2 are ordered such that Im > |^2| > Im {^s} |^4|- 

Notice that Im 2| > 0, whereas Im |fe3,4| < 0. At (3 = 0, the wavenumbers are independent 
of 0. As /? increases from zero, the dependencies of the wavenumbers upon 9 are observed to 
be highly asymmetric with respect to 9 = Tr/2. In the limit [3 ^ 1, the ^-dependencies of 
the real parts of the wavenumbers become antisymmetric relative to ^ = 7r/2, whereas the 9— 
dependencies of the imaginary parts of the wavenumbers become symmetric relative to = 7r/2. 

In relation to planewave propagation, a topic of considerable current interest is whether the 
phase velocity is negative or positive (Lakhtakia et al. 2003). Negative phase velocity (NPV) 

is closely related to the phenomenon of negative refraction (Ramakrishna 2005). Planewave 
propagation with NPV in S is signified by (Mackay & Lakhtakia 2004b) 



Re{^} •P(r,a;) < 0, 



(48) 



where P(r, w) is the time-averaged Poynting vector; conversely, positive phase velocity (PPV) 
in S is signified by 

Re{k} •P{r,uj) > 0. (49) 
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In S', NPV is signified by Re {fe'} • P'(r',u;') < and PPV by Re {M} ' P!{r',io') > 0. Issues 
concerning NPV propagation for isotropic cfiiral mediums (Mackay 2005) and FCMs arising as 
homogenized composite mediums (Mackay & Lakhtakia 2004b) have been reported previously. 

For the medium of interest here, NPV propagation occurs in E' provided that (Mackay 2005) 
Re|v^eV)MV) -?'('^')} x 



Re I ^/€'{LJ')^J,'{uJ') + ^'{io') } X Re 

For the same medium, by virtue of (48), NPV propagation occurs in S provided that 

n{k)<0, (51) 

the form of the real-valued NPV parameter Cl(k) being provided in Appendix 1. 

Let us return to the numerical example considered in Figures 1 and 2. The /3^-regimes of NPV 
and PPV, as determined by evaluating fl{k) for k G ^2, ^3, ^4|, are mapped in Figure 3 
with respect to the relative speed /3 G [0,1) and wavevector orientation angle 9 G [0, tt). The 

medium clearly docs not support NPV propagation when /3 = 0; i.e., all plane waves in E' must 
be of the PPV kind. As P increases, the /3^-regimes supporting NPV propagation in S emerge 
in the range 7r/2 < ^ < tt for the relative wavenumbers ki^2, and in the range < ^ < 7r/2 for 
the relative wavenumbers ^3^4. 

Finally in this section, we note that an alternative derivation of the dispersion relations (42) 
in S may be developed via the Lorentz transformation of the corresponding dispersion relations 
in S'. NPV arises when this Lorentz transformation brings about a change of sign in the angular 
frequency. 




(50) 



4 Dyadic Green functions 

The problem of finding the (frequency-domain) field phasors generated by a specified distri- 
bution of sources within a linear medium is conveniently tackled by means of dyadic Green 
functions (DGFs) (Tai 1994). For the isotropic chiral medium in reference frame S', the DGFs 
are well-known (Lakhtakia 1994b). However, explicit representations of DGFs are generally 
unavailable for anisotropic and bianisotropic mediums (Mackay & Lakhtakia 2006). In this 
section we exploit the constitutive relations derived in §2, together with the Beltrami-like fields 
introduced in §3, to establish a convenient spectral representation of the DGFs for the FCM 
described by (29). 

Let a source electric current density phasor J_g{r,u) and a source magnetic current density 
phasor J„j(ii, w) exist, from the perspective of the non-co-moving observer. Extending the 
approach adopted in §3 wherein Beltrami-like fields are introduced to aid the planewave analysis 
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in S, we recast Jemfcj^) ^ Beltrami current density phasors 



Wi{r,u) = ^ [iv'{oo')J,{r,u;)-J^{r,u;) 



(52) 



The Beltrami-like field phasors generated by the source terms (52) may then be expressed in 
terms of the DGFs G (r — s,(jj') as 



Q/l, u^) = Q^fc, ^) + J^g^ir- s, to') • W^(r, co) d^s, (^=1,2), 



(53) 



where V is the region containing the source current density phasors. The complementary func- 
tions Q 2 ^) given by 



CO 



(54) 



wherein £^(r , u>) and ^(r, u>) satisfy the relations 

V X E{r, id) - iuoB{r, w) = 1 

V X H_{r, uj) + iuDir, w) = J ' 

along with 

Dir, u) = €o g{uj') ' E{r, to) + i,/e^^{io') • H_{r, uj) 
Ml, w) = {{uj') • E{r, u) + Ho fJ.{uj') • H{r, lu) 

The DGFs in (53) are provided as the solutions of the differential equations 

V X (r - s, J) + {-if ko UJ) 'GAr- s, J) = 8{r - s) /, (£ = 1, 2) 



(55) 



(56) 



(57) 



with 5( • ) being the Dirac delta function. By implementing the spatial Fourier transforms 

(9> = l^t fc' (-^^ • ^) ' = 1' 2) (58) 

with (57), the components of the spectral DGFs G" {q,uj') with respect to the Cartesian basis 
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emerge as 

g\{q,u; 

g,iq,uJ 

&{iq,LO 
g^{q,LU' 

g\{q,^' 
g\{q,cv 



11 



12 



13 



21 



22 



23 



31 



32 



33 



{q ' x)'^ — i^iiii 



Ai 



(q ' yf - K^Kf 



[{q 'x){q'y)+ ihl {hf + g' z)] Ai 

[q ' x[q ' z_ + K^) — iR\q • y] Ai 

[{q 't){q'y}- if^l + i)] 

Ai ) , 

[q'y{q'z + Kl) + iK{q • x\ Ai 
[q't{g'z + kl) + ik\q • y] Ai 
\_q ' y {q ' z_ + Ki) — ik\q • ^] Ai 

+ k{ + q'z)[k\-k{-q' z)] Ai 



+ {k{f + 2k{klq' z - {k\f , 



£2(^' 

£2(^''^ 
£2(«''^ 



11 



12 



13 



21 



22 



23 



31 



32 



33 



A, 



[- {q 'x){q'y)+ ikl {-k^ + g * l)] A2 

\q'x{q'z-kl)+ ik\q • y] A2 

[- {q 'x) {q'y) - ikl {-kl + g • 1)] A2 



- [q-yf + i^\~K.2 

\q'y_{S_'z-kf)- ik^q ' x] A2 
\q'x{q'z-k^)- ik\q • y] A2 
\q'y_{q'z-kl)+ ik\q • A2 
[(k*2 - + g • i) (4 + - g • z)] A2 



A, 



(g • 1) - (g • x) - (g • J - K*2) - g • g {kl + 
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Having established the spectral DGFs G"^ ^ (^j'^Oj obtain the S field phasors generated by 
the source phasors Jgmfc'^) ^ 



■i 



g[{q,u')e^p[iq' {r-s)] dq\ 



j ds 

(63) 



and 

1 



47^3 



V \\Jq 



{q, w') exp [ig' (r - s)] dq 



■^7^ |^£''i (^,w')exp [ig • (r-s)] dgj • [W {oj') JeiL,^^) - ImiL:'^)]^ ds. 

(64) 



5 Discussion 

The Tellegen constitutive relations for an isotropic chiral medium moving at constant velocity 
are presented in (23) and (29). The availability of these constitutive relations facilitates a full 
analysis of the planewave characteristics of the medium, and also enables the spectral DGFs to 
be derived in a convenient form. In contrast to earlier studies, the analysis presented herein is not 
restricted to low relative speeds (Hilhon 1993; Ben-Shimol & Censor 1995, 1997). Furthermore, 
the constitutive relations (23) and (29) establish that the uniformly moving isotropic chiral 
medium in fact belongs to the category of FCMs. 

In §3, the analysis of planewave propagation in reference frames S' and S, is aided by the 
introduction of the Beltrami field phasors Q' {jdi^') in (31) and the Beltrami-like field phasors 
Qi 2 fc''^) i^ (36)) respectively. They facilitate a decouphng of the Maxwell curl postulates. A 
key property of Q'^ ^ (i'j is that the curl of Q'-^ ^ (r', lo') is a scalar multiple of Q'^ ^ (r', w'), as 
demonstrated in (34). Such fields are known as Beltrami fields and their properties are firmly 
established (Hillion k, Lakhtakia 1993; Lakhtakia 1994b, b). In contrast, the curl of Q^^ {l^^) 
is not generally parallel to ^ iz^^), as may be observed from (38). The Beltrami-like field 
phasors Q , „ (r, u) therefore represent an important extension of the usual Beltrami field concept 
which can be traced back to at least the late 1880s (Beltrami 1889; Silberstein 1907; Trkal 1919). 

Owing to their relatively large parameter space, linear bianisotropic mediums support a richer 
palette of planewave properties than do anisotropic and isotropic mediums, as has been high- 
lighted lately by investigations of NPV propagation (Mackay k, Lakhtakia 2004b) and optical 
singularities (Berry 2005). The planewave study presented in §3 reveals that the bianisotropic 



14 



FCM described by the constitutive relations (29) generally supports four independent wavenum- 
bcrs for each direction of propagation, from the perspective of a non-co-moving observer (the 
exception being propagation perpendicular to the direction of translation for which only two 
independent wavenumbers are supported). This constrasts with the two independent wavenum- 
bers supported by the isotropic chiral medium from the perspective of the co-moving observer. 
We see in Figure 3 that an isotropic chiral medium, which does not support NPV propagation 
from the perspective of a co-moving observer, does support NPV propagation from the per- 
spective of a certain class of non-co-moving observers, provided that the relative speed is suffi- 
ciently high. This finding is consistent with results presented for an isotropic dielectric-magnetic 
medium moving at constant velocity (Mackay & Lakhtakia 2004a). This is also consistent with 
a study which showed that a FCM arising as a homogenized composite medium can support 
NPV propagation provided that the gyrotropy parameter of the gyrotropic constituent medium 
is sufficiently large (Mackay Sz Lakhtakia 2004b). 

While explicit representations of DGFs are available for isotropic mediums, these are generally 
not available for anisotropic and bianisotropic mediums (Mackay & Lakhtakia 2006). However, 
as shown in §4, the field phasors for the FCM described by the constitutive relations (29) may 
be formulated in terms of spectral DGFs. By exploiting the constitutive relations (29) and 
the Beltrami-like field phasors Q 2 ^ convenient representation of the spectral DGFs is 

estabUshed in (58)-(62). 
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Appendix 1 



As indicated by (51), whether or not plane waves propagate with negative phase velocity in a 
Faraday chiral medium is determined by the sign of the parameter ^}{k). The form of Q{k) is 
provided here; a complete description of the derivation of ri(fc) is available elsewhere (Mackay 
& Lakhtakia 2004a). 

As a function of the relative wavenumber k, 0(fc) may be expressed in terms of the constitutive 
parameters (25)-(27) and the wavevector orientation angle 9 as 



Re 



X Re 



/x^*(w') 



sm( 



[/x**(c.')]'-[m^*(^')]' 



k* sine- IT* 
k* I [i/*{io') (lap + 1) + ifj,<'*{uj') {a - a*)] cos^ 9 



+\Tff/*ico')sm^9- [iJ*{uj'){a*T + aT*)+iiJ,3*ico')ir-T*)] sin^cos^ 
+ [h**{lo')C^*{u') - hS*{lo')C**{uj')] [{\af + l) cos0 - a*T sm9] 
-i [/x**(a;')C**(a;') - ijP*{J)i3*{u')\ [{a -a*)cos9-T sin 9] 



(65) 
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Herein, 



with 



a 



-^13 ~ -^11-^33 
L12L23 — LisL 



22 



Lu = e\u;') + 



^22 



L13L23 + L12L33 



) 1 5^ — 



i^33 = e'{io') 



L12 = 



Ll3 = 



A; sin^ 



-^23 = i 



k sin 9 



(66) 



and r = ^9{Lo') + kcose. 
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Figure 1: Real (left) and imaginary (right) parts of the constitutive parameters %*, ^ind 
X^, where % G {e,^, /u}, plotted against the relative speed /? = v/cq. Solid curves represent Xj 
dashed curves represent Xzi ^'^d broken dashed curves represent Xg- Constitutive parameters 
in S': e' = 6.5 + ^' = 1 + iO.2, and n' = 3.0 + iO.5. 
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Figure 2: Real (left) and imaginary (right) parts of the four relative wavenumbers ^1,2,3,4 plotted 
against the relative speed /3 = f/co and the propagation angle 6 = cos~^^ • k (in degree). 
Constitutive parameters in S': e' = 6.5 + il.5, ^' = 1 + iO.2, and fi' = 3.0 + iO.5. 
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Figure 2: continued 
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Figure 3: /30~regimes of NPV and PPV, as viewed from the inertial frame S, with respect to the 
relative speed P = v/cq and the propagation angle 9 = cos~^ v • k (in degree), as determined 

by the parameter i^(k) for the four relative wavenumbers k E ^2) ^3) ^4|- Constitutive 

parameters in E': e' = 6.5 + il.5, ^' = 1 + iO.2, and /x' = 3.0 + iO.5. 
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